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Abstract 

We show that flavor evolution of a system of neutrinos with continuous energy spectra as in su- 
pernovae can be understood in terms of the response of individual neutrino flavor-isospins (NFIS's) 
to the mean field. In the case of a system initially consisting of fg and i>e with the same energy 
spectrum but different number densities, the mean field is very well approximated by the total 
angular momentum of a neutrino gyroscope. Assuming that NFIS evolution is independent of the 
initial neutrino emission angle, the so-called single-angle approximation, we find that the evolution 
is governed by two types of resonances driven by precession and nutation of the gyroscope, respec- 
tively. The net fiavor transformation crucially depends on the adiabaticity of evolution through 
these resonances. We show that the results for the system of two initial neutrino species can be 
extended to a system of four species with the initial number densities of Ve and i>e significantly 
larger than those of I'x and v'x- Further, we find that when the dependence on the initial neutrino 
emission angle is taken into account in the multi-angle approximation, nutation of the mean field 
is quickly damped out and can be neglected. In contrast, precession-driven resonances still govern 
the evolution of NFIS's with different energy and emission angles just as in the single-angle ap- 
proximation. Our pedagogical and analytic study of collective neutrino oscillations in supernovae 
provides some insights into these seemingly complicated yet fascinating phenomena. 
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I. INTRODUCTION 



Experiments on solar, atmospheric, reactor, and accelerator neutrinos have found over- 
whelming evidence for neutrino flavor transformation via vacuum oscillations and the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect (see [l| and references therein). The parame- 
ters governing these phenomena have largely been determined except for the neutrino mass 
hierarchy, the vacuum mixing angle 613, and the CP- violation phase. For supernova neu- 
trinos, in addition to the MSW effect induced by neutrino forward scattering on electrons 
in matter 2, sl, novel phenomena of flavor transformation can occur (see 1^ and references 
therein) due to neutrino forward scattering on other neutrinos IJH]. This is because near 
the proto-neutron star produced in a supernova, the number density of neutrinos can exceed 
that of electrons. The dominance of neutrino-neutrino over neutrino-electron forward scat- 
tering couples flavor evolution of neutrinos that have different energy and travel on different 
trajectories, thereby giving rise to new phenomena. Assuming no CP-violation in the neu- 
trino sector, js-lO] first showed that collective oscillations of supernova neutrinos can occur 
for the vacuum-mass-squared difference Sml^ relevant for atmospheric neutrino oscillations. 
This collective flavor transformation can be treated effectively as 2 x 2 mixing between z/g 
and (z/g and Ux), where (i^x) is a superposition of and Ur {v^ and Vr)-, because z/^ 
and Vr i^^i and Ur) have identical emission characteristics and interactions in supernovae. 



Subsequent to the study in [10| , a number of authors have investigated collective oscilla- 
tions of supernova neutrinos (see ^ for a review). The most prominent feature of collective 
oscillations is the stepwise swap of the Uf, and z/^. energy spectra at some characteristic "split" 



energy 
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12| . As this spectral swap for supernova neutrinos depends on the neutrino mass 
hierarchy (i.e., the sign of Sml^) and 613, neutrino signals from a Galactic supernova can 
potentially provide a sensitive probe of both these unknown parameters. 

Recently, it was found that spectral swaps for supernova neutrinos can occur at more than 



one energy 13l-ll5|. Specifically, using supernova neutrino emission parameters with initial 
number densities of Ue and z/g significantly larger than those of Ux and Ux, jisl] found that 
for an inverted mass hierarchy (IH, the lightest mass eigenstate being dominantly Ux), an 
additional swap occurs in the spectra of Ue and Ux at a rather low energy. In contrast, using 
very different parameters with initial number densities of z/^, and Ux significantly larger than 
those of z/g and Ue, [l4| and [l5| showed that spectral swaps occur at more than one energy 
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for both IH and normal mass hierarchy (NH, the heaviest mass eigenstate being dominantly 

In this paper we focus on the occurrence of spectral swaps for the case where the initial 
number densities of Uf. and Ue are significantly larger than those of and u^. As a concrete 
example, we assume that neutrinos are emitted in one of the flavor eigenstates Ua (repre- 
senting any one of Ue, u^, z/g, and Ux) from a sharp neutrino sphere of radius -R^, = 10 km 
with a normalized spectrum of the form 

1 98 F''^ 

USE) = ^7Tr^exp(-4E/(E.J), (1) 

where {Ey^) is the average energy. We take {E^^) = 10 MeV, {EpJ = 15 MeV, and {E^J = 
(Ep^) = 24 MeV. We further take the neutrino luminosity to be L^^ = 10^^ erg s~^ for each 
species so that the initial number densities of z/g, Ue, i^x, and are in ratios of 2.4 : 1.6 : 1 : 1. 

We define = I^Jri^gl > and 6^ to be the parameters for vacuum mixing between 
z/g and ipe and Px), where < 6^ < 7r/4 or 7r/4 < 6^ < n /2 corresponds to the NH or 
IH, respectively. In accordance with the experimental limits on ^^13, we take ^ 1 (NH) 
or 6'v = 7r/2 — ^ 1 (IH). Following the definition of |l^, we represent each neutrino by 
a neutrino flavor isospin (NFIS) s of magnitude 1/2 in a three-dimensional flavor space. In 
particular, 

et/2, for z/e or z/^,, 
-e[/2, for Vx or z?e, 

where is the unit vector in the 2;-direction of the flavor space. It is convenient to use 

{5m? /2E, for z/p and z/^., 
/ , 
—5m?'/2E^ for u^, and Ux-, 

to label the NFIS of an initial z/^ with energy E' as s^. 

We follow the flavor evolution of neutrinos in the above example from the neutrino sphere 
to a radius r = 250 km, where oscillations have effectively ceased. We focus on the IH with 
6m? = 3 X 10~^ eV^ and assume 9^ = 10~^ in view of matter suppression of mixing, but 
otherwise ignore the matter above the neutrino sphere. We further assume that the flavor 
evolution of a z/^ emitted in a non-radial direction is identical to that of another Ua with the 
same energy but emitted radially. This is the so-called "single-angle" approximation. The 
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FIG. 1: An example of stepwise spectral swaps for the IH and for the case where the initial number 
densities of Vf. and Ve are significantly larger than those of and Vx- The neutrino flux spectra 
at the neutrino sphere (dashed curves) and at r = 250 km (solid curves) are shown. The dashed 
curves with larger maxima are for the initial Ve and Ve while those with lower maxima are for the 
initial and Vx- It can be seen that the and Vx spectra are swapped for E ^ 7.6 MeV and the 
Ve and Vx spectra are swapped for £^ > 3.1 MeV. 

effective total number density at radius r > Ry for neutrinos emitted initially as Va is taken 
to be 

The evolution of the NFIS with time t, or equivalently, with radius r is governed by 



d ^ 
~dt 



d 

( 

dr 



Jo 



(5) 



where = — sin2^ve^ + cos26've^ = — sin26've^. — cos2^ve^, Gp is the Fermi constant, 
and u' runs over all NFIS's as summation over a and integration over E are performed. For 
convenience, we adopt the natural units, where the Planck constant h and the speed of light 
c are set to unity, in writing expressions and equations throughout the paper. 
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FIG. 2: The swap factor fs{uj,t) at r = 250 km as a function of 2uj/6m'^ for the example shown 
in Fig. [TJ Note that the negative range of the horizontal axis corresponds to and Ux while the 
positive range corresponds to Ve and Vx- 

Figure [T] compares the neutrino flux spectra at the neutrino sphere (dashed curves) and 
at r = 250 km (solid curves). The vertical scale measures the relative energy-differential 
fluxes, which are oc L^^f^^{E)/ {E^J at the neutrino sphere. It can be seen that the 
and Ux spectra are swapped for E > 7.6 MeV and the Ue and Ux spectra are swapped for 

^ 3.1 MeV. To highlight these spectral changes, we define a "swap factor" 

fsiu;,t)^sUt)/si,M, (6) 

where sl^ ^(t) = S(^(t) • e^. The value of fs{uj,t) = 1 or — 1 corresponds to complete survival 
or conversion of the initial Ua, respectively. The swap factor fs{u},t) at r = 250 km for the 
above example is shown as a function of 2u/5m?' in Fig. [2j 

The main goal of this paper is to explain how spectral swaps occur in the above example. 
Our approach is pedagogical. In Sec. [Ill we consider the net effect of all the NFIS's in the 
above example as a "mean field" to which an individual s^^ responds. We show that the 
evolution of this mean field can be well described by the motion of a "neutrino gyroscope" 
over a wide range of neutrino densities. In Sec. IIIII we discuss in detail the precession and 
nutation of the neutrino gyroscope. In Sec. IIVI we show that the response of an NFIS to 
the mean field is characterized by two types of resonances driven by the precession and 
nutation of the neutrino gyroscope, respectively. These resonances and the adiabaticity of 
evolution through them give rise to the stepwise spectral swaps. In Sec. |Vl we relax the 
single-angle approximation and discuss the evolution of NFIS's along different trajectories in 
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realistic supernova environments. We show that the trajectory-dependent evolution strongly 
suppresses nutation of the neutrino gyroscope. Consequently, the spectral swaps in this case 
are determined by precession- driven resonances only. We summarize our results and give 
conclusions in Sec. IVll 



II. MEAN FIELD OF NFIS'S AND NEUTRINO GYROSCOPE 

The NFIS for an initial z/g and that for an initial with the same energy (or uj) are equal 
in magnitude but opposite in direction [see Eq. ([2])]. It can be seen from Eq. ^ that these 
two NFIS's remain equal in magnitude but opposite in direction during their subsequent 
evolution. The same is also true of the NFIS's for an initial z/g and an initial with the 
same energy. Hereafter, s^i for w > or w < refers to the NFIS for an initial z/g or z/g, 
respectively. Then Eq. (j5]) can be rewritten as 

d 



dr 



/oo 
g{uj')s^iduj' 
-co 

where /i(r) = 2\/2G pny^ir) . In the above equation. 



(7) 



UAE^)-c.J,SE^), fora;>0, 
9{u.) = ^ X <! (8) 
CefuAEuj) - Csfp^iE^), for u <0, 



where = n^,{R^)/n^^{R^) = {E^J/{E,yJ and E^ = 6m^/2\uj\. 
We define 

/oo 
9ico)s^du (9) 
-oo 

and consider it as the mean field representing the net effect of all NFIS's to which an 
individual s^^ responds. Clearly, an exact description of the mean field requires solving 
Eq. ([7]) for all NFIS's. However, we can give an approximate description based on the 
evolution of a small number of NFIS's. To see this, we show 2g{uj) / 5im? as a function of 



2uj / 5im? for our supernova example in Fig. |3l There are three "spectral crossings" IJ] at 
(jj = w_, 0, and respectively, for which g{u) = 0. In each of the four spectral regions 
separated by these crossings, the magnitude of g{u)) is large only for a relatively narrow 
range of u. 

We consider that the four spectral regions (— oo, (w., 0), (0, co+), and oo) can be 
represented by four effective NFIS's Sg°\ s^^"*, si°\ and si^"*, respectively, which correspond 
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FIG. 3: Example supernova neutrino spectrum shown in terms of 2g{ijj)/5m? as a function of 
2u!/5m?'. There are three spectral crossings (filled circles) at u; = a;_, 0, and a;+, respectively, for 
which g(ijj) = 0. 



to a t^e, a z/a;, a i^a;, and a z/g at the neutrino sphere. We then approximate S by 

where the evolution of each si'^'* is governed by 

d 



dr 



[a;,H.-/x(r)S(°)]. 



In the above equations, 



a. 



1 



g{u)\du, 
uj\g{uj)\duj, 



(10) 



(11) 



(12) 
(13) 



and other quantities are defined similarly. For our supernova example, ae = 0.78, ax = 0.20, 
as = 0.10, and a^ = 0.35, while 2uja/5m^ = 0.15, 0.032, -0.027, and -0.11 MeV^^ for si°\ 



si^"*, s^^"*, and s^^', respectively. 

At the neutrino sphere, S = S'-^-' = As can be shown from Eqs. ([7]) and ffTTl) . the 

components of S and S*^°^ parallel to are conserved during the subsequent evolution. We 
numerically obtain their components perpendicular to H^, S± (solid curve) and 5*^"'' (dashed 
curve), and show these as functions of jj,{r)/uje in Fig.HJ It can be seen that S"^^ closely tracks 
S± for fi{r)/u)e ^ 20 but large deviations occur for fi{r)/u)e < 10. In particular, S^^^ diverges 
from S± at fj,{r)/ue < 4. The above results can be understood by comparing /x(r)|S| with 
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FIG. 4: Comparison of S± (solid curve) and S]_ (dashed curve) as functions of fi{r) /oOf, for the 
supernova example. 

the spread in co for the spectrum shown in Fig. [3l For fi{r)/cOe ^ 20, /i(r) |S| ^ 3ue exceeds 
the spread over the entire spectrum, which is ~ 2ue- So all NFIS's evolve collectively in this 
regime and the four effective NFIS's included in 8'-°^ are sufficient to give a good description 
of S. For fi{r)/ue ~ 10, /i(r)|S| ~ Ue approaches Aue = 0.66a;e, which is the spread in the 
spectral region (a;+, oo) calculated from 



Consequently, for fi{r)/ue < 10, the NFIS's in the above spectral region are no longer well 
represented by Se''^ and large differences between S^^^ and S± occur. Eventually, none of the 
effective NFIS's can represent their respective spectral regions and 5'^°'' diverges from S±. 



Eq. for the evolution of s^, and use the results to obtain a better approximation for S from 
Eq. (Q. This procedure may be repeated until successive approximations for S converge. 
While this approach does not save numerical efforts compared with solving Eq. ([7]) directly, 
it motivates an analytic study based on the zeroth order mean field S'-^-*, especially when 



paper. We first consider a simpler case and then apply the results from this case to discuss 
the supernova example in Sec. |Vl 




(14) 



In a formal approach, we can use S*-"-* as the zeroth order approximation for S to solve 



S^^-* can be understood with simple models. We carry out such a study in the rest of the 
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A. System Initially Consisting of v^, and Vf, Only 



In our supernova example, the initial number densities of z/g and Ue are significantly larger 
than those of and z/^. To facilitate an analytic study, we consider a simpler system initially 
consisting of v^. and only. We take normalized emission spectra of the form in Eq. ([T]) 
with {E^^) = {EpJ = 12 MeV and luminosities L^^ = 1.2 x 10^^ and Lp^ = 0.8 x 10^^ erg s"^ 
so that the initial number densities of Ue and z/g are the same as those in the supernova 
example. In this case, the effective neutrino spectrum g{uj) reduces to 



6m^ \fuAEuj), forw>0, 
2^ 



M-)~x^---' ' (15) 



afp^{E^), for c; < 0, 

where a = ripXRu) / nyXRy) = 2/3. 

We consider that the zeroth order approximation for the mean field S = h(co)Si^duj 
is given by 

S(o)^sf)+asf , (16) 
where s^^'' and S2*''* are the NFIS's for an initial z/g and an initial z/g with 

POO 

ui = ujh{uj)duj, (17a) 



JO 

1 



UJ2 = — ujh{u)du, (17b) 
a J~oc 

respectively. Specifically, ui = —uj2 = /iv and 2^1^ /dm? = 1/9 MeV""*^. We show S±_ (solid 
curve) and sf^ (dashed curve) as functions of ^{r)/^^ in Fig. [51 It can be seen that the 
comparison between S±_ and 5*^"'' is very similar to that in the supernova example [note that 
the horizontal scales for Figs. H] and O are related by /i(r)//iv = 1.35/i(r)/u;e]. 

The swap factor fs{uj,t) at r = 250 km for the system initially consisting of z/g and z/g 
only is shown as the solid curve in Fig. [61 The swap factor f^g\u,t) calculated from 

|;s^ = s^x [a;H.-/x(r)SW], (18) 

which uses S^^-* to approximate S, is shown as the dashed curve. It can be seen that just 
as in the supernova example, fs{uj,t) and ff\uj,t) have two characteristic split energies, 
one in the region of w > and the other in the region of w < 0. Ahhough fPico,t) has 
a different split energy from that of fs{uj,t) in the region of < 0, they have the same 
quahtative behavior, especially the same spht energy in the region of > 0. Our goal is to 
understand the behavior of fg^\uj,t) analytically. 
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FIG. 5: Comparison of S± (solid curve) and (dashed curve) as functions of ii{r)/ii^ for the 
system initiahy consisting of Ve and Ve only. 
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FIG. 6: Comparison of the swap factors fs{uj,t) (solid curve) and fg'\io,t) (dashed curve) at 
r = 250 km for the system initially consisting of and z^e only. 



B. Neutrino Gyroscope as the Approximate Mean Field 

For the system initially consisting of Uf. and z/g only, the evolution of s^'''' and Sg''"' is 
governed by 



|-s(°) = srx [-;.3v-Mr)S(°)]. 



(19a) 
(19b) 



For convenience, we will drop the superscript "(0)" but otherwise use the symbols in the 
same meaning as in the above equations. It is useful to consider the time evolution of Si 
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and S2 at a constant fi governed by 



di 

si = —Si = Si X (/ivHv - /iS), (20a) 
S2 = ^S2 = S2 X (-/ivHv - /iS). (20b) 



As discussed in 
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17l | and repeated below, the system governed by the above equations is 
mathematically equivalent to a gyroscope in a uniform gravitational field. 
From Eqs. f l20ap and ( ]20b|) it is straightforward to show that 

Q = /iS X Q, (21a) 

S = X H„ (21b) 

where 

Q = si - as2 + — H,. (22) 

From Eq. ( 12 la^ it can be shown that Q = |Q| is conserved. With the definition of a unit 
vector f = Q/Q, Eqs. (12 lap and ( ]21b[) can be rewritten as 

S = -f X £• + af , (23a) 

S = -f X g, (23b) 
A'' 

where a = S-r and g = fifi^QH^. The above equations describe the motion of a gyroscope in 
a uniform gravitational field with an acceleration of gravity g (see Fig. Uh)- The gyroscope 
has a spinning point particle of mass 1/^ attached to the tip of a massless rod of unit 
length. The spin of the gyroscope is along the direction f of the rod and it can be shown 
from Eqs. (12 lap and ( I21bp that the magnitude a of the spin is conserved. As can be seen 
from Eq. ( ]23bp . the component of the total angular momentum S parallel to Hv is also 



conserved. 

If fi varies smoothly with time, the neutrino gyroscope described above evolves through 
a series of configurations corresponding to a continuous range of fi. The motion of this 
gyroscope provides a well-studied mechanical analog to the evolution of the approximate 
mean field. 
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FIG. 7: Illustration of (a) the neutrino gyroscope and (b) the Euler angles and <j) used to describe 
its motion in Frame I. The gyroscope has a spinning point particle of mass 1/^ attached to the tip 
of a massless rod of unit length. The spin is along the direction f of the rod and has a magnitude 
cr. The gyroscope is in a uniform gravitational field with an acceleration of gravity g. The wiggly 
curve on the spherical surface crudely indicates the trajectory of the tip of f as the gyroscope 
executes precession (change in cp, arrow labeled by (pp) and nutation (change in 6) at a constant fx. 

III. PRECESSION AND NUTATION OF THE NEUTRINO GYROSCOPE 

We specify the unit vector f of the neutrino gyroscope by two of the Euler angles, 6 and 
0, in Frame I: 

f = sin 6' sin (pe], — sin 9 cos + cos 9 e\, (24) 

where ej., e^, and are the unit vectors associated with the three axes of Frame I and 
= — Hv (see Fig. [7]). We focus on the IH with 9^ <^ 1, for which = —sm29^el. — 
cos26've^ ^ —el- Thus, ^ is the upward direction and g (in the same direction as Hy) 
points downward (see Fig. [7^). When the neutrino gyroscope starts at yU ^ /i^, Q ~ Si — as2 
points in the direction of e^, i.e., it is in the upright position. This initial configuration can 
give rise to interesting subsequent evolution. 
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A. Motion of the Neutrino Gyroscope at Constant /x 



We first discuss motion of tfie neutrino gyroscope at constant /i. Using Eq. f l23ap . we can 
express S in terms of 9 and as 

Sx = — cos (j) — (j)sm.6 cos 9 sin 0) + cr sin 9 sin 0, (25a) 
/i 

Sy = — (6^ sin + sin 9 cos 9 cos 0) — cr sin 6* cos 0, (25b) 
/i 

5*2 = —0 sin^ 6' + cr cos 6* = — cos 26'v, (25c) 

wliere we liave used conservation of Sz in the exact form (i.e., no approximation made for 
^ 1) in the last equation. In addition to a and S^, the third conserved quantity of the 
gyroscope is its total energy: 

^gyro = - if • g, (26a) 

= ^(^9^ + 02 sin^ + + fi^Q cos 9. (26b) 

Conservation of -Egym can be shown using Eqs. f l23ap and ( ]23b[) . We can also derive an 
explicit equation of motion from Eq. ( I23bl) : 

9 — 0(0 cos 9 — fia) sin 9 = fifi^Q sin 9. (27) 

In general, both 9 and of a gyroscope evolve with time, and the corresponding motion 
(see Fig. [7]) is referred to as nutation (9) and precession (0). For the neutrino gyroscope, 
conservation of a, Sz, -Egyro can be combined to give a conserved effective energy associated 
with nutation only: 



where 



Eg = E-gyro - f = ^ + ^efr(^), (28) 

/^x li (Sz — (7 COS 9Y ^ „ / X 

2 sm 9 

As shown in Sec. IIIIBt nutation of the neutrino gyroscope mostly occurs around the min- 
imum of the effective potential Ves{9). This potential minimum corresponds to 9 = 9p, for 
which 



/i 



d9 

eS 



d9^ 



{Sz - acos9p){a - SzCOs9p) n ■ a n fin \ 

ji ^ /ivtt/ smfc'p = 0, (30a) 

sin 9p 

^p sin^ 9p + ^20p cos 9p- fia^ =ujI>0. (30b) 
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In Eq. (]30bll . (/)p is the instantaneous precession frequency at the potential minimum and 
can be obtained from Eq. fl25cp after Op is solved from Eq. (I30al) . Note that Eq. (I30al) is the 
same as 

COS 9.p - /iO") + /iyUvQ = 0, (31) 

which can also be obtained by setting = and = 0p in Eq. ( 1271) . In other words, the 
minimum V^s corresponds to the maximum 9 [due to conservation of Eq, see Eq. (l28l) ]. and 
hence ^ = 0. 

To the leading order, nutation can be approximated as oscillations of 9 around 9p in 
response to the potential 

V,M-V,^{9p) + '^^{9-9pf. (32) 
In this approximation, the evolution of 9 can be described by 



9 ^ 9p- r]cos{ujnt + /3), 
9 !=a TjUJn sin(co'„t + /3), 



(33a) 
(33b) 



where rj is the amplitude of oscillation (or nutation) and /3 is a constant phase. The values 
of rj and /3 are determined by the initial values of 9 and 6* at t = 0. Using Eq. (I25cl) to relate 
at 6* to 0p at 9p^ we obtain to the first order in 77, 



ixa — 20p cos 9p 



sin 9n 



cos(w„t + /3), 



(j) ^ + (j)pt - 7] 



fia — 2(j)p cos 9p 
LOn sin 9p 



sin(w„t + (3), 



(34a) 
(34b) 



where 7 is a constant determined by the initial value of at t = 0. 
From Eqs. fl25all -fl25cD and flMijl - flMbl) we obtain 



Sx + iSy 



— COS 9jj — a \ sin 9„ + - — 
/i 



2 fi 



^ _ 20P - fiS^ I 



+ 1 



bp - jj,Sz 



-i(Unt+l3) 



OJr, 



(35) 



which is accurate to the first order in rj. The above equation shows that to the leading order, 
the angular momentum S can be described by three precessing vectors with different ampli- 
tudes and different precession frequencies, which are 0p, 0p + cOn, and (pp — cOn, respectively. 
This leading-order expression of S provides a very useful analytic description of the motion 
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of the neutrino gyroscope. Without loss of generahty, hereafter we set the phase constants 
/3 = 7 = in Eq. §5^. 



B. Motion of the Neutrino Gyroscope for Slowly Decreasing fi 

We now extend the discussion in Sec. IIII Al for a constant /i to the case where ^(t) slowly 
decreases with time from a large initial value /i(0) ^ fi^. Specifically, we consider the system 
that initially consists of monoenergetic u^. and z/g only and is governed by Eqs. fl2Uap and 
f l20b|) . For this system, 2fi^/6m^ = 1/9 MeV-\ a = 2/3, and 

-\ 2 



= /i(0) 




r 



t + T 



(36) 



where /i(0) = 2.52 x 10^/iv and r = 8.45//iv. The above form of /i(t) corresponds to n^^{r) in 
the supernova example. Taking 6^ = 10^^ (IH), we numerically solve Eqs. f l20ap and (120bp . 
Using the instantaneous Si and S2 from the numerical results, we construct a gyroscope at 
each specific value of /x. We show cos 6* (dashed curve) for the series of gyroscopes as a 
function of /i(t)//iv in Fig. [HI Using a and Q for each gyroscope and the fact that Sz is 
conserved even when fi changes with time [see Eq. fl21b|) ]. we construct an "instantaneous" 



Vcs{0) from Eq. (!29|) . We define ^max and 6'inin as the solutions to Ves{0) = Eg, where Eg is 
the corresponding instantaneous effective energy [see Eq. (l28i) ]. We show the evolution of 
cos 6'max and cos as the dotted curves in Fig. [HI We also calculate 9p corresponding to 
the minimum of the instantaneous Ves{0) and show the evolution of cos 6p as the dot-dashed 
curve in Fig. [8] Finally, we calculate (pp and Un from Eqs. fl25cl) and fl30bp and show 4>p/ f^v 
(dashed curve) and w„//iv (sohd curve) as functions of nit)/ jj,^ in Fig. [9] 

It can be seen from Fig. [8] that the general trend of 9 follows the evolution of dp. In other 
words, nutation of the neutrino gyroscope occurs around the minimum of the instantaneous 
Vcs{d) as /i(t) slowly decreases from a large initial value yu(0) /iv It is also clear that 
the motion of the gyroscope falls into two distinct regimes separated by a critical /Xcr ~ 
119/iv For n > ficr, the nutation amplitude is extremely small and 6* ~ 6'p to very good 
approximation. As fi drops below fi^r, & initially stays small even as 9p increases. This can 
be understood from Fig. [9], which shows that the nutation frequency a;„ is small at yU ~ ^„ 
and is practically zero at Ha- As ji decreases further, ujn becomes sufficiently large and 6 
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FIG. 8: Nutation of the neutrino gyroscope shown in terms of cos as a function of fi{t) / fi^ (dashed 
curve) for = 10~^ (IH). The gyroscope is characterized by 2n^/5m? = 1/9 MeV-\ a = 2/3, 
and = /i(0) |l - Y^l - [T/(t + r)]2|^, where /i(0)//Uv = 2.52 x 10^ and r = 8.45//^^. The 
dot-dashed curve is for cos^p, which corresponds to the minimum of the instantaneous 14fj(^), and 
the dotted curves are for 

cos ^rnax snd cos^min, which Correspond to V^f[(^) — Eq. 




200 150 100 50 
l^(t)/Hv 



FIG. 9: Dimensionless precession frequency (t)p/^iv (dashed curve) and nutation frequency uJn/ t^-v 
(soHd curve) as functions of fi{t)/fi^ for the neutrino gyroscope shown in Fig. [8j 

starts to oscillate around 6p. The amplitude of this oscillation is also that of nutation and 
can be taken as r] ^ (^max — ^mm)/2. The longer Un stays small at /i ~ the larger i] is 

for /i < /icr- 

We have constructed the series of gyroscopes for specific values of fi using the instan- 
taneous Si and S2 numerically calculated for the /i(t) in Eq. ( 136|) . In fact, so long as fj,{t) 
slowly decreases from some large initial value /i(0) ^ /iv, the characteristics of such gyro- 
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scopes essentially depend on the values of /i but not the specific functional form of /i(t). To 
see this, we consider the set of parameters a, Q, and 6p that characterize the gyroscope at a 
specific value of /i. We can choose three equations to solve for these parameters as follows. 
From the definitions of a, S, and Q, we obtain 

= - (37a) 

S' + Q^ + 2^Q cose+(^y = (37b) 



Applying Eqs. f l25cp and fl37bp to ^ = 6p, at which 6 reaches its maximum value of (^)r 



TjoUn [see Eqs. (I33ap and (]33b|) ]. we further obtain 



— sin^ 9p + a cos 6p = S^, (38a) 

+ — sm^9p + a +Q +2 — Qcos6'„+ — ^ — - — . (38b) 

\ fi J \ fi fi \fi J 2 



Treating 77 as a small parameter and ignoring the 77^ term in Eq. (]38bp . we can solve this 



equation along with Eqs. fl37ap and fl38ap to obtain a, Q, and 6p [note that (pp is given in terms 
of cr, Q, and 6p by Eq. (1311) : see Appendix [X] for a different but equivalent method to obtain 
these parameters]. As this approximate solution assumes 17 = 0, it is a "pure-precession" 
solution for which the gyroscope always stays at the minimum of the instantaneous 

For 6'v ^ 1, the pure-precession solution gives 

Op < 1, (39a) 
1 + a 
1 — a 

1 — a / 1 + a\ , , 

-/i - 2 fi,, (39c) 



'p 7/iv, (39b) 



2 " \l-a, 

at yU ^ (see Appendix [B] for more detailed discussion of the initial motion of the neutrino 
gyroscope), and 

Op vr, (40a) 
0P ^ yUv, (40b) 
Un — 2/iv, (40c) 
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at /i ^ In general, the 9p, (j)p, and Un calculated for the pure- precession solution are 
within ~ 1% of the values shown in Figs. [8] and [91 

Equation fl39cl) suggests that Un becomes very small as decreases to some critical value 
/icr- In Appendix [Cl we show that 

(1 - ^/a) 

which is /icr ~ 119/iv for a = 2/3, in excellent agreement with Figs. |8]and|9l At /i = ficr, we 
have (see Appendix ICj) 

^"'^'^ ii + ^yi^^'' ' ^^^^^ 

' (42b) 

(42c) 

2 1 — \/« 



As ^ 1 at ^ and at /i = /^cr for ^ 1, the neutrino gyroscope stays in the upright 
position and behaves like a sleeping top (e.g., [igI at /i > /^cr (see Fig. [S]). 



IV. RESONANCES DRIVEN BY THE NEUTRINO GYROSCOPE 

In this section we return to the system initially consisting of and z/g with spectra of the 
form in Eq. ([1]). As discussed in Sees. Ill Al and III Bt the approximate mean field of NFIS's 
for this system can be described by the neutrino gyroscope. We now try to understand 
the evolution of an individual NFIS in the system in terms of its response to the neutrino 
gyroscope. Specifically, we study the evolution of governed by 

^s, = X [wH, - /i(t)S] = [uje\ + /i(t)S] X s,, (43) 

where S is the total angular momentum of the neutrino gyroscope discussed in Sec. IIIIBI 



A. Precession-Driven Resonance 



We first ignore nutation and consider only precession of the neutrino gyroscope. With 
= (and 7 = as noted in Sec. IIII Al) . Eq. ( 135]) becomes 

Sx + iSy ~ iS±_ exp{i(j)pt), (44) 
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where 

S± = \ — cos 6p — a I sin6'p. (45) 




Equation fl44j) represents a vector rotating in the xy-plane of Frame I. Let Frame II rotate 
with an angular velocity (pp^]. relative to Frame I. Then S is a non-rotating vector in Frame 
II and can be chosen as 

S ^ S^el + S^e];, (46) 

where e^^ is the unit vector in the y-direction of Frame II (see Fig. [TU]) . We rewrite Eq. f l43p 
in this frame as 

= [{uj + fiS, - <j)p)el + fiS±ef] x = Hn x s^. (47) 

Note that in the above equation u and Sz are constants but 0p and S± are functions of 
/i(t). A resonance occurs when the z-component of Hn vanishes. We refer to this as the 
precession-driven resonance. We denote the value of ^ at which an individual NFIS s^^ goes 
through this resonance as ^res,pi which can be obtained from 

UJ = 0p(/Xres,p) — l^res,pSz- (48) 

Here 0p(/ires,p) is the value of 0p at /i = /ircs,p- Using the neutrino gyroscope in Fig. [HI we 
show /ires,p//^v as a function of w/yUv (solid curve) in Fig. [H] (note that the pure-precession 
solution gives essentially the same result). 

The evolution of in Frame II based on Eq. fHTj) is very similar to the MSW effect 
(see Fig. [T2|) . At a specific time t, s^^ precesses around the instantaneous Hn with an 
angular velocity Hn. For slowly varying /i(t), the evolution of is adiabatic in that the 
precession adjusts to the instantaneous angular velocity (the direction and magnitude of 
which are both changing slowly in general) but the angle between s^^ and Hn remains fixed. 
Therefore, the initial z/g or z/g represented by remains in the same flavor following adiabatic 
evolution if the initial and final directions of Hn are the same, but is fully converted into 
a or if the initial and final directions of Hn are opposite. For the neutrino gyroscope 
under consideration, Hn ~ fiSzel at t = corresponding to fi ^ /iv, where we have used 
fiSz ^ 0p, fJ'\Sj_\ in this limit. At large times corresponding to /i <^ /iv, Hn ~ (w — /iv)^^, 
where we have used <pp = fi^ in this limit. Consequently, the initial and final directions of 
Hn are the same for u > fi^ but are opposite for u < fi^. In the latter case, a precession- 
driven resonance occurs when the ^-component of Hn vanishes before the direction of Hn is 
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Frame I Frames II & III 




FIG. 10: Illustration of relations among Frames I, II, and III. Frame I is fixed to the laboratory. 
Frame II has the same z-axis as Frame I and rotates around this axis with an angular velocity 
(j)pe\ relative to Frame I. When only precession of the neutrino gyroscope is considered, the net 
effective field interacting with an individual NFIS has two components in Frame I: one fixed 
in the direction of and the other rotating with an angular velocity (/>pe^,. It is convenient to 
use Frame II to discuss the precession-driven resonance as the net effective field becomes a non- 
rotating vector Hii in this frame. Frame III has the same x-axis as Frame II but its z-axis is in 
the direction of Hjj, which makes an angle 6 with respect to e^. Frame III is used in discussing 
the nutation-driven resonance. See text for details. 




FIG. 11: Values of Hres,p/ fJ-v (solid curve) and /ires.n/^iv (dashed curve) as functions of lo/hv An 
NFIS 

^ui goes through the precession-driven (nutation-driven) resonance at // — /^res,p (A'resjn) s-s 
the neutrino gyroscope in Fig. [8] evolves through /i(t). 
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FIG. 12: Illustration of nearly full conversion of an initial Ve represented by with uj < fi^ fol- 
lowing adiabatic evolution through a precession-driven resonance. In Frame II, adiabatic evolution 
corresponds to precession of around Hn with a fixed angle between the two vectors. At /U ^ 
Hn nearly coincides with e\. At /_i = ^res,p corresponding to the resonance, the z-component of 
Hii vanishes. At ^ /^v, Hn is opposite to e\ for uj < Consequently, the initial and final 
directions of Hn are nearly opposite, and the same is true of s^^. 

reversed (see Fig. [T2|) . Thus, when only precession-driven resonance matters and adiabatic 
evolution applies, an initial z/g with u > fi^ remains as a i/g, while an initial z/g with u < fi^ or 
an initial Ue (with w < 0) is fully converted into a or Ux, respectively. This is basically the 



explanation for the stepwise spectral swap originally discovered in Uj (see also discussion 
in Q). 

For adiabatic evolution, the rate at which the direction of Hn changes must be slow 
compared with the precession frequency of s^: 

d H„ 



dt I Hii I 



r^"""^"" « |H„|. (49) 
nil 



The above condition is most stringent at resonance when the 2;-component of Hn vanishes 
and I Hii I becomes very small (see Fig. [T2i) . We define the adiabaticity parameter for this 
precession-driven resonance as 

_ I Hii X dUu/dtlres _ \d{^p - iJ,S,)/dt\ 

' = Wis " ii^' • ^ ^ 

Adiabatic evolution obtains for Ap ^ 1 (note that Xp is defined differently from the usual 
adiabaticity parameter for the conventional MSW effect). Using the neutrino gyroscope in 
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FIG. 13: Adiabaticity parameters Xp (solid curve) and A„ (dashed curve) for evolution through 
precession-driven and nutation-driven resonances, respectively, as functions of /i(t)//iv for the neu- 
trino gyroscope shown in Fig. [HI 

Fig. [HI we show Xp as a function of ^{t)/fi^ (solid curve) in Fig. [131 It can be seen from this 
figure that evolution through the precession-driven resonance is adiabatic at /i < llSyUv, but 
is extremely non-adiabatic at /i ^ ~ 119/iv. As /i ^ fia- corresponds to the sleeping-top 
regime with ^ 1, the small values of S± oc sin6'p [see Eq. ( I45l) ] result in Ap ^ 1 in this 
regime. 

B. Nutation-Driven Resonance 

Now we consider both precession and nutation of the neutrino gyroscope. The terms 
proportional to r] in Eq. (!35l) contain the factors exp(iu;„t) and exp{—iuJnt), which corre- 
spond to rotation with angular velocities UnG\ and — t<j,„e^, respectively, relative to Frame 
II. However, frames rotating with these angular velocities are not convenient to use because 
fi{t)S±e^y, and hence Hn, rotate in such frames. To find the appropriate frames, we first 
consider Frame III with its axes defined by the unit vectors (see Fig. [TOl) 





cos 6e„ — sin 6e 



(51b) 



(51a) 
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where 



cos 6 



sin 6 



iHiil 



IHtt 



(52a) 
(52b) 



Note that just hke Frame II, Frame III also rotates with an angular velocity (ppe], relative to 
Frame I. Using Eq. (1351) (with /3 = 7 = as noted in Sec. IIII A[) . we write S in Frame III as 



S ^ SA + S^e" + r]^ sin S cos(w„t)e™ 

/i 



+ 77- 



bp2(j)p fiSz . . .N^iii , 0p f / 

- sm(co'„t)e^ + 77— cos b cos(w„t)e . 



(53) 



The last two terms in the above expression can be rewritten as two vectors rotating with 
angular velocities w„e™ and — u;„e^^\ respectively, relative to Frame III: 



{S^ + iS')iu = i 



v4^ 

2 /i 



cos 6 — 



(54) 



As we will see shortly, a new resonance occurs for I Hill = Using Eqs. ( l30b|) . (EI]), 
f l38ap . and (HTll . we can rewrite the above resonance condition as: 



{u + 0p)(w + 2^Sz - 30p) = 0. 



(55) 



The term with the factor exp(— 2a;„t) in Eq. fl5^ vanishes for u = —(pp, while that with the 
factor exp{iunt) vanishes for u = 3(f)p—2fiSz. We will see that the new resonance corresponds 
to u = —(pp. So we can ignore the term with the factor exp{—iunt) in Eq. ( 15^ when treating 
this resonance. We choose Frame IV to rotate with an angular velocity Wne™ relative to 
Frame III. The term with the factor exp(iu;„t) in Eq. (15^ represents a vector parallel to the 



unit vector e^^ in the ?/-direction of Frame IV. In this frame Eq. fH3|) effectively becomes 



dt 



-H-iil — Un + 1](Pp Sm cos Unt 1 6^ H I COS I e. 



UJr, 



X s„ 



(56) 



It can be seen that a new resonance indeed occurs when |Hii| = a;„ if we ignore the small 
contribution proportional to r] in the term associated with in the above equation. We 
refer to this as the nut at ion- driven resonance because it is driven by the nut at ion- dependent 
component of S. We denote the value of /i at which an individual NFIS s^^ goes through 



this resonance as /ires.n, which can be obtained from 



(pp(^f^res,n) • 



(57) 
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Here (j)p{fires,n) is the value of (j)p ai fi = /Xres.n- Using the neutrino gyroscope in Fig. |H1 
we show /ires,n//^v as a function of u/^^ (dashed curve) in Fig. [11] (note again that the 
pure-precession solution gives essentially the same result). 

To see that u = 30p — 2fiSz, which also gives |Hn| = w„, does not correspond to a 
resonance, we recall that the term with the factor exp{iunt) in Eq. ( 15^ vanishes for this u. 
We choose Frame V to rotate with an angular velocity — a;„e^^^ relative to Frame III and 
rewrite Eq. ( l43i) in Frame V effectively as 

r]4>p I 20P - fiS^ 



d 



|^|Hn| +Un + sm () cos w„t J + I cosd H I 



(5^ 



where is the unit vector in the y-direction of Frame V. It can be seen that the term 
associated with in the above equation never vanishes, and consequently, there is no 
resonance for u = 34>p — 2fiSz. 

The adiabaticity for evolution through the nutation-driven resonance can be discussed 
similarly to the case of precession-driven resonance studied in Sec. IIV A[ The evolution of 

in Frame IV is governed by 

d 

— = Hiv X s^. (59) 

At resonance Hiv is given by 

\ J 

We define the adiabaticity parameter for the nutation-driven resonance as 

_ |Hiv X c/HivMlres \d{\Hii\ - Un) / dt\res /^,x 

" = PfTls ~ ITT 12 • voij 

I Hiv I res l-H-Ivlres 

In Eqs. fl60|) and fl6T|) . we have neglected oscillatory terms proportional to r]. Using the 
neutrino gyroscope in Fig. El we show A„ as a function of yu//iv (dashed curve) in Fig. [131 It 
can be seen from this figure that evolution through the nutation-driven resonance is adiabatic 
at /i//iv ~ 10-110 and becomes non-adiabatic outside this range. In particular, evolution 
is extremely non-adiabatic at /i ^ /icr ~ 119/iv and at /i ^ 5/iv as the very small nutation 
amplitude r] (see Fig. [8]) results in A„ 1 in these two regimes. 



C. Evolution through Resonances Driven by the Neutrino Gyroscope 



Based on the discussion in Sees. IIV Al and IIVB| an NFIS may experience two types of 
resonances driven by precession and nutation of the neutrino gyroscope, respectively. For the 
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NFIS S(^, a precession-driven resonance occurs at a; = (j)p{fires,p) — fJ'res,pSz, and a nutation- 
driven resonance occurs at w = — 0p(/ircs,n)- Formally these two resonances coincide at 
/^res,p = yUres,n = /^cr for u = = -0p,cr ~ (1 + A/a)/iv/(l " a/o) [see Eq. fl42bp]. Noting 



that 0p — )■ (1 + — a) at /i ^ /icr [see Eq. (]39b|) ] and 0p — > /iv at /i ^ /iv [see 



Eq. ( ]40b[) ]. we introduce loa = —(1 + a)/Uv/(l — a), ujb = — /Wv, and uc = /Uv to define ranges 



of u with different resonances. It turns out that for u = ud slightly larger than ujc-, there 
are two possible values for /Ures,p- Altogether, the varieties of resonances experienced by s^^ 
can be classified into six categories: 

I. for uj < Ucr, experiences only a precession- driven resonance at /ires.p > A^cr! 

II. for oJcT < (jJ < uja, experiences a nutation-driven resonance at fires,n > A^cr, then 
a precession-driven resonance at /ires,p < /^cr, and finally a second nutation-driven 
resonance at /i^^^ „ < fi^es,p] 

III. loT ua < < ujb, experiences a precession-driven resonance at /ires,p < /^cr followed 
by a nutation-driven resonance at fires,n < /^res.p; 

IV. for cus < a; < cuc s^^ experiences only a precession-driven resonance at ^res,p < yUcr! 

V. for < < wd, experiences two precession-driven resonances at /ires,p ^ /^cr and 
/^res,p < /^res,p, respectively. 

VI. for a; > ud, s^^ does not experience any resonance. 

For the neutrino gyroscope in Fig. [HI ficr ~ 119/iv, i^cr ~ — lO/iv, ~ — 5/iv, and w^:) ~ 
1.5/Xv. The six categories of resonances for this example are shown in the top left panel of 
Fig.[Il 

Only precession-driven resonances are involved for u in ranges I, IV, and V, and there are 
no resonances for u in range VI. We first consider the overall evolution of s^^ for u in these 
ranges using the neutrino gyroscope in Fig. [81 As nutation is unimportant for these cases, we 
focus on Hn as the net effective field interacting with s^^ [see Eq. ( HT]) ]. At /i = /i(0) ^ 
Hii is essentially in the direction of e], and is either aligned (initial u^, u > 0) ot anti- 
aligned (initial z/e, < 0) with Hn. We show the subsequent evolution of relative to 
Hn by comparing cos6'^ = s^^ ■ e^js^l (dashed curve) with cos6'n = Hn • ey|Hn| (solid 
curve) in Fig. [TH For u/fi^ = —12 in range I, is initially anti- aligned with Hn (cos^^^ ~ 
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FIG. 14: Ranges I- VI of uj based on occurrences of precession-driven (solid curve) and nutation- 
driven (dashed curve) resonances (top left panel) and example evolution of cos 6^) (dashed curves) 
and cos^ii (solid curves) as functions of ^{t)/fi^ for cj//iv = —12 (I), —0.5 (IV-a), 0.5 (IV-b), 1.1 
(V), and 3 (VI), respectively. See text for details. 



— 1 but cos^ii ~ 1). At /i//iv ~ 122, cos^n vanishes and a resonance occurs. However, 
evolution through this resonance is extremely non-adiabatic (see Fig. [T3|) . Consequently, 
cos6i_j is unaffected while cos6'ii changes drastically from ^ 1 to ~ — 1 immediately after the 
resonance. Subsequent evolution of s^^ is essentially adiabatic with cos 9^ oscillating around 
cos ^11 and eventually settling to ~ —1 again. This kind of evolution applies to all u in 
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range I, for which there is no net flavor transformation. For w//iv = —0.5 (0.5) in range IV, 
there is a resonance at yu/yUv ~ 36 (24). Evolution through the resonance is adiabatic (see 
Fig. [T5]) and s^^ stays anti-ahgned (ahgned) with Hn during the entire evolution. As the 
initial and final directions of Hn are opposite, there is full fiavor conversion for uj in range 
IV. For u;//iv = 1-1 in range V, there are two resonances at /i//iv ~ 14 and 1, respectively. 
Evolution through both resonances is essentially adiabatic (see Fig. [13]) and stays aligned 
with Hn during the entire evolution. As cos^^n changes sign twice, the initial and final 
directions of Hn are the same and there is no net fiavor transformation for uj in range V. 
Finally, for a;//iv = 3 in range VI, there is no resonance and evolution is adiabatic. So cos 9^^ 
oscillates around cos^^n? indicating that is always aligned with Hn. There is no net fiavor 
transformation for uj in range VI. 

Resonances driven by both precession and nutation of the neutrino gyroscope are involved 
for UJ in ranges II and III. We discuss the evolution of s^^ for these ranges using Hjv as the 
net effective field. We define cos6'iv = Hiv • e!,/|Hiv|. Neglecting terms proportional to 
Tj, we obtain cos6'iv ~ sgn(|Hn| — Ci;„)cos6'n (see Sec. IIVBI) . where sgn(x) is the sign of 
X. Using the neutrino gyroscope in Fig. [HI we compare the evolution of cos^^^^ and cos6'iv 
for u in ranges II and III in Fig. [151 For bjj^i^ = —6 in range II, is initially aligned 
with Hiv. A nutation-driven resonance occurs at /i//iv ~ 210. However, evolution through 
this resonance is extremely non-adiabatic (see Fig. [T3|) . So cos^^^^^ is unaffected although 
cos^^iv jumps from ^ — 1 to ^ 1 immediately after the resonance. Then cos6'iv vanishes at 
/i//iv ~ 99 corresponding to = 30p — 2iiSz and |Hn| = Wn, but this is not a resonance (see 
Sec. lIVBl) . A precession-driven resonance occurs at /i//iv ~ 86 and a second nutation-driven 
resonance occurs at /i//iv ~ 43. Evolution through both these resonances is adiabatic (see 
Fig. [13]). Consequently, at /i//iv < 99, cos6'tj oscillate around cos6'iv, eventually settling to 
~ —1 again. There is no net flavor transformation for bj in range II. For u/ = —4 in range 
III, cos^iv vanishes at /i//iv ~ 89, which is not a resonance. A precession-driven resonance 
occurs at at /u//iv ~ 69 and a nutation-driven resonance occurs at yu/yUv ~ 19. Evolution 
is adiabatic throughout and there is no net flavor transformation. For uo / = —2 also in 
range III, the evolution of cos ^^iv is similar to that for ojj = —4. However, for ojj = —2, 
evolution through the precession-driven resonance at /u//Xv ~ 51 is adiabatic while that 
through the nutation-driven resonance at /i//iv ~ 4 is non-adiabatic (see Fig. [T3]) . This kind 
of evolution applies to —3 < u;//iv < ~1 and results in signiflcant flavor transformation. 
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FIG. 15: Ranges II and III of w based on occurrences of precession-driven (solid curve) and 
nutation-driven (dashed curve) resonances (top left panel) and example evolution of cos 9^1 (dashed 
curves) and cos^iv (solid curves) as functions of ^{t)/^^ for uj/^^ = —6 (II), —4 (Ill-a), and —2 
(Ill-b), respectively. The dotted curve in the top left panel shows the value of ^/^v for which 
u) = 3(j)p — 2^Sz as a function of oj/fiy. Note that no resonance occurs at this although the 
corresponding cos 6yv vanishes. See text for details. 

In summary, evolution of can be understood in terms of the resonances it experiences. 
A resonance does not affect the net flavor transformation if evolution through it is very 
non-adiabatic. Net full flavor conversion results from adiabatic evolution through an odd 
number of resonances while little net flavor transformation results from adiabatic evolution 
through an even number (including zero) of resonances. As discussed in Sees. II V A] and HV Bl 
evolution through a resonance driven by either precession or nutation at yU > /icr is extremely 
non-adiabatic. In contrast, evolution through a precession-driven resonance at /i < /icr is 
essentially always adiabatic. We introduce a parameter yUtr to discuss the adiabaticity of 
evolution through a nutation-driven resonance at /i < ^cr'- the evolution is adiabatic (non- 
adiabatic) when such a resonance occurs at yUtr < yU < yUcr (/^ < /^tr)- We choose /itr to 
correspond to an adiabaticity parameter A„ = 0.5. For the neutrino gyroscope in Fig. |8l 
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/itr ~ 12/iv and with u ~ — 3/iv goes through a nut at ion- driven resonance at /i = /itr (see 
Figs. H] and [13]). 

Now the swap factor shown as the dashed curve in Fig. E] can be understood based 
on the above discussion and the occurrences of resonances hsted in the beginning of this 
subsection and shown in Figs. [H] and US This curve is re-plotted as the sohd curve in 
Fig. [T6b with ranges I to VI for u indicated. Recall that ff^ ^ — 1 corresponds to nearly 
full flavor transformation. This applies to u in range IV, for which there is only a single 
precession-driven resonance at /i < and evolution through this resonance is adiabatic. 
A variety of evolution results in f^^^ ^ 1 corresponding to little net flavor transformation 
for Lo in ranges I, II, V, and VI: non-adiabatic evolution through a single precession-driven 
resonance at /i > /icr (1)? non-adiabatic evolution through a nutation-driven resonance at 
yU > /icr followed by adiabatic evolution through a precession-driven resonance and a second 
nutation-driven resonance aX n < /icr (H), adiabatic evolution through two precession-driven 
resonances at /i < /icr (V), and adiabatic evolution with no resonance (VI). A precession- 
driven resonance and a nutation-driven resonance occur at /i < /icr for uj in range III and 
evolution through the precession-driven resonance is always adiabatic. However, evolution 
through the nutation-driven resonance is adiabatic only for u < —3/iv in this range and 
becomes more and more non-adiabatic as u increases above ~ —3/iv. Consequently, a 
transition from f^^ ^ 1 towards —1 occurs at a; ~ —3/iv in range III. 

D. Application to System of Neutrinos with Continuous Spectra 

As discussed in Sec. Ill Bt the total angular momentum of the neutrino gyroscope, which 
we denote as S^^^ again for clarity, approximates the mean field S of the NFIS's in the system 
of neutrinos with continuous spectra as specified in Sec. Ill A[ The oscillations of S"^^ and 
S± shown in Fig. [5] reflect the nutation of the gyroscope. It can be seen from this figure that 
large deviations of S'^°^ from S occur only at /i//iv < 20, where the nutation amplitude of S 
rapidly decreases. In addition, S± sharply drops at /i//iv < 2. The small nutation amplitude 
of S at /i//iv < 20 affects the adiabaticity of evolution through the nutation-driven resonance 
for u in range III (see top left panels of Figs. [Ml and ITSl) . In fact, the evolution is extremely 
non-adiabatic for —3 < u/fi^ < —1. On the other hand, evolution through the precession- 
driven resonance is adiabatic for these values of u, which results in net full flavor conversion. 
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FIG. 16: Comparison of the swap factors fs{uj,t) (a) and f^\LO,t) (b) at r = 250 km for the 
system initially consisting of and Ue only as in Fig. [6l The solid, dashed, dotted, and dot-dashed 
curves are for 6^ = 10~^, 10^^, 10^^, and 10^^, respectively. Ranges I- VI of w shown in (b) are 
based on occurrences of precession-driven and nutation-driven resonances (calculated for 9^ = 10^^ 
but vahd for (9^ < 1). 

Thus, compared with the results based on S''^^ the swap factor is ~ —1 for a wider range 
of oj as shown by the solid curve in Fig. [T6k . In principle, the sharp decrease of S± at 
/i//iv < 2 could affect the adiabaticity of evolution through the precession- driven resonance 
at the lower /i for u in range V (see Fig. [T4|) . However, in practice this has little effect 
(see the solid curve in Fig. [T6k) as resonances at such low values of fi only occur for a very 
narrow range of u and adiabaticity is affected for an even narrower range of u. 

To further illustrate how adiabaticity of evolution through precession-driven and nutation- 
driven resonances affect net flavor transformation, we increase 6^ from 10~^ to 10^'^, 10~^, 
and 10~^, respectively. For a larger 9^, S± in the sleeping-top regime of /i > is larger as 
the initial 6p of the neutrino gyroscope becomes larger [see Eq. fH5l) and Appendix[B]. On the 
other hand, the nutation amplitude rj is smaller at /i < ficr as it grows less at ;U ~ due to 
a shorter nutation period ~ 2ti /ujn^cr for a larger 6^ [see Eq. fl42cl) ]. Consequently, evolution 
through a precession-driven resonance at yU > /icr becomes less non-adiabatic while that 
through a nutation-driven resonance at /i < /Xcr becomes more non-adiabatic. The former 
effect becomes quite large for 9^ = 10^^ as partial flavor conversion occurs for u in range I 
(dotted curves in Fig. [16]), while the latter effect is already significant for 9^ = 10^^ (dashed 
curves in Fig. [T6l) as more flavor transformation occurs for u in ranges II and III relative to 
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the case of 9^ = 10~^ (solid curves in Fig. [T6ll . For 6^ = 10~^, S± is sufficiently large initially 
and 7] remains small at all /i. Consequently, evolution through precession- driven resonances 
is adiabatic for u in ranges I to V while nutation-driven resonances have no effect on the net 
flavor transformation. This can be seen from the dot-dashed curves in Fig. [161 which show 
that net full flavor conversion occurs for u in ranges I to IV with a single precession-driven 
resonance but there is no net flavor transformation for u in ranges V and VI with two and 
zero precession-driven resonances, respectively. 

V. COLLECTIVE NEUTRINO OSCILLATIONS IN SUPERNOVAE 

In this section we consider the system of neutrinos exhibiting the spectral swaps shown in 
Figs. [Hand [21 As described in Sec. [H this system initially consists of u^, i^e, ^x, and with 
continuous spectra. In addition, the initial number densities of z/g and are significantly 
larger than those of and Ux- We ffist show that the swap factor shown in Fig. [2 can be 
understood in terms of the six different kinds of flavor evolution that have been discussed 
in Sees. IIV CI and II V D I for the system initially consisting of z/g and z/g only. The evolution of 
cos9^ corresponding to Fig. [2lis shown in Fig. [T7]for 2uj/5rn? = —0.90, —0.50, —0.32, 0.02, 
0.14, and 0.50 MeV~^, respectively. It can be seen that these six kinds of evolution are very 
similar to those shown in Figs. [TH and [151 for the six ranges of u discussed in Sees. IIV CI As 
the comparison of S^^-* and S shown in Fig. [4] for the system of four initial neutrino species is 
similar to that shown in Fig. [5lfor the system of two initial neutrino species, the differences 
between the evolution based on S*^''-' and S are also similar to those discussed in Sec. IIVDI 
Therefore, we conclude that the flavor evolution of a system with initial number densities 
of and Ue significantly larger than those of and can be understood in terms of the 
resonances driven by precession and nutation of a neutrino gyroscope. 

Next we consider the flavor evolution of the system exhibiting the spectral swaps in 
Figs, m and [21 by relaxing the single-angle approximation used to produce these results. In 
the so-called "multi-angle" approximation, neutrinos are emitted from the neutrino sphere 
with equal probability in the forward directions, which are defined to be < O^m < 7r/2. 
Here 6'em is the angle with respect to the radial direction at the point of emission. Under the 
multi-angle approximation, an NFIS can be specified by the corresponding neutrino energy 
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FIG. 17: Example evolution of cos^o; as a function of ^{r)/(jje for the system shown in Figs. [T] 
and El for 2uj/5m'^ = -0.90 (I), -0.50 (II), -0.32 (III), 0.02 (IV), 0.14 (V), and 0.50 (VI) MeV-^ 
respectively. Note the similarity to the evolution shown in Figs. [TH and [TJ] for the six ranges of u 
based on occurrences of precession-driven and nutation-driven resonances. 

and emission angle in terms of u and e = cos6'em- The evolution of s^^^ is governed by 

?2 /"oo /•! 



d 

dr 



-U^-2fx{K)^ / gico')dco' / 

' J-oo Jo 



1 



- 1 



(62) 
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where fi{Ru) = '2'\/2Gpn^^{R^), n^^{R^) is given by Eq. (jl]) for r = R^, and 



DJr) 



1 - (1 - e^)Rl/r^. 



(63) 



Equation (!62|) reduces to Eq. ([7]) under the single-angle approximation that S(^ ,;(r) = 

Sa;,e=l(?^)- 

Using the same neutrino emission parameters as for Figs. [H and El we follow the flavor 
evolution of the system with initial number densities of z/g and Ue significantly larger than 
those of Ux and Ux under the multi-angle approximation. The angle-averaged swap factor 
(/s((X', e, r))^ at r = = 250 km is shown as a function of 2u/5w? in Fig. [T8k . This average 
factor is defined as 

jlfs{uj,e,T)[e/D,{r)]de _ fs{u,e,r)[e/D,{r)]de 



{fs{uj,e,r)). 



Io{^/De{r)]de 



'l-Rl/r^ r^/Rl 



(64) 



and can be used to calculate the effect of flavor transformation on e.g., neutrino reaction 
rates at radius r. Compared with the swap factor shown in Fig. |2] for the single-angle 
approximation, {fs{co,e,rf))^ is identical for w > but shows large deviations for u < 0. 
On the other hand, swap factors /s(w, e, rj) for specific values of e shown in Fig. [T8b have 
the same general structure as shown in Fig. HJ This can be understood from the precession- 
driven resonance as discussed below. 



It is convenient to define 



g{uj)duj 
g{uj)duj 



^iO^€^d€^ 



DJr) 



de. 



(65a) 
(65b) 



and rewrite Eq. ( 1621) as 

d 
dr 



X 



H, 



2/i(i?.)^' 



DJr] 



(66) 



D^{r) ' ' ' 

Our numerical results show that both Sa and S;, precess around Hv with the same frequency 
0p(r) at a specific r for R^ < r rf 



This is consistent with the conclusions of [18 1 
based on symmetry arguments. The synchronized oscillations of Sa,x/Sa,± and Sb,x/Sb,± due 
to precession are shown for r = 60-80 and 200-250 km in Figs. [T9k and b, respectively. 
Here the x-axis is in the plane perpendicular to and the subscript "±" denotes the net 
perpendicular component. It can be shown from Eq. (|66|) that 5*^,2 = 8^-6^ (e^ = — H^) 
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FIG. 18: Angle-averaged swap factor {fs{uj,e,rf))^ (a) and angle-specific swap factors fs{u},e,rf) 
(b) at Vf = 250 km for the system with the same neutrino emission parameters as for Fig. [2] but 
in multi-angle approximation. From left to right, the curves for a; < in (b) are for e = 10~^, 
0.24, 0.49, 0.74, and 0.99, respectively. The sharp changes in /^(a;, e, rj) occur at 2u/Sm^ = -1.63, 
-1.54, -1.24, -0.86, -0.24, and 0.14 MeV-\ respectively, for which fs{uj,e,rf) = (filled circles). 



is conserved (see Fig. [20b). Components of the "mean field" Sb/D^{r) — Sa are shown as 
functions of r for e = (solid curve) and 1 (dashed curve) in Figs. [20b and c. Clearly, the 
mean field experienced by ^ is different for different e. It can also be seen from Fig. [2Ub 
that the perpendicular component Sb,±/D^{r) — Sa,± increases sharply at r^r ~ 64 km just 
as S± does at the critical point in the single-angle approximation (see Fig. [4|). However, in 
contrast to the single-angle approximation, nutation of the mean field is damped out very 
quickly in the multi-angle approximation, and therefore, can be neglected. 

Based on the above discussion, we only need to consider the precession- driven resonance 
in explaining the swap factor fs{u}, e, Vf). A precession-driven resonance occurs when the z- 
component of the net effective field interacting with ^ vanishes in the co-precessing frame 
(see Sec. IIV Ap . It can be seen from Eq. (!66|) that s^^^ goes through a precession- driven 



resonance when 



2 



r 



DJr 



(67) 



However, evolution through a resonance at r < is extremely non-adiabatic and does 
not result in any net flavor transformation. Thus, we expect that net full conversion 
[/s(w,e,r/) = -1] occurs only for Wres(e,rcr) < u < ujres{e,rf). At r = r^r, a;res(e, ^^cr) < 
for < e < 1 and is larger for a larger e (more radial trajectory, see Fig. [20H). At r = Vf, 
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FIG. 19: Synchronized oscillations of Sa,x/Sa,± and Sh^x/Sb,± due to precession for the system 
shown in Fig. [18] for r = 60-80 (a) and 200-250 km (b). In either panel two curves are shown but 
they are indistinguishable. 



Co 



0.086 


." ' 1 1 


"'1 1 

(a); 




0.008 p 


0.085 


1 

:\ 




N 

1 


0.006 I 


0.084 






u 

N 

si 
CO 


0.004 - 
0.002 if 


0.083 


..1 1 1 1 ' 




^ 



50 100 150 200 250 
r (km) 




100 150 
r (km) 



250 



to 



-D 
CO 




> 



6 
to 



S 




100 150 200 250 
r (km) 



FIG. 20: Components of the mean field as functions of r: (a) Sa,z (solid curve) and Sb^z (dashed 
curve), (b) Sb^z/ D^{r) — Sa,z for e = (solid curve) and 1 (dashed curve), and (c) Sb^±/ D^{r) — Sa,± 
for e = (solid curve) and 1 (dashed curve), (d) The value of 2(^i-es(e, r^r) / StvP at r^r ~ 64 km shown 
as a function of e. From left to right, the filled circles correspond to 2(^i-gs(6, r^r 

)/(5m2 -1.57, 

-1.50, -1.24, -0.81, and -0.21 MeV"\ respectively. These values are very close to those values 
of ^ijjjbrr? that correspond to the filled circles at cj < in Fig. 118b . 
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D^{r) ^ 1 and Sh,z/ D^{r) — Sa,z becomes very small (see Fig. [20b). so Wresle, ff) ~ (i>p{rf) and 
is nearly independent of e. Using our numerical results for (ppir), Sa,z, and Sb^z, we obtain 
2wres(e,rcr)/5m2 ^ -1.57, -1.50, -1.24, -0.81, and -0.21 MeV^^ for e = lO^^, 0.24, 0.49, 
0.74, and 0.99, respectively, and 2u)j-cs{i,rf) /Sm^ ^ 0.14 MeV~^ for < e < 1. These results 
are in excellent agreement with Figure [T8b . 



VI. CONCLUSIONS 



Using a system initially consisting of z/g and Ue with the same energy spectrum but 
different number densities, we have shown that flavor evolution of this system in the single- 
angle approximation can be understood in terms of the response of individual NFIS's to the 
mean field, which is very well approximated by the total angular momentum of a neutrino 
gyroscope. The evolution of an NFIS is governed by two types of resonances driven by 
precession and nutation of the gyroscope, respectively. A resonance does not affect the net 
flavor transformation if evolution through it is extremely non-adiabatic. Nearly full flavor 
conversion occurs following adiabatic evolution through an odd number of resonances but 
there is no net flavor transformation following adiabatic evolution through an even number 
(including zero) of resonances. The detailed results on NFIS evolution are presented in 
Figs. [H [151 and dSl and discussed in Sees. [IVC] and [IVDl 

We have also shown that the above results for the system of two initial neutrino species 
can be extended to a system of four species with the initial number densities of and z/g 
significantly larger than those of and Ux- Further, we find that when the multi-angle 
approximation is adopted instead of the single-angle approximation, nutation of the mean 
field is quickly damped out and can be neglected. In contrast, precession- driven resonances 
still govern the evolution of NFIS's with different energy and emission angles just as in the 
single-angle approximation. These results are presented and discussed in Sec. [Vl 

In conclusion, we have presented a detailed analysis of collective neutrino oscillations 



in supernovae for the case where the initial number densities of anc 



larger than 



studies 



20 



;hose of and u^. We note that some earlier works (e.g.. 



are significantly 



19|) and two recent 



21| have similar goals to ours but used very different methods. Our approach 



is mostly pedagogical and analytic. It is our hope that along with other parallel efforts, we 
have provided some insights into the seemingly complicated yet fascinating phenomena of 
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collective neutrino oscillations. 
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Appendix A: Parameters of the Neutrino Gyroscope at a specific fx 

In addition to the procedure given in Sec. IIII Bt the parameters of the neutrino gyroscope 
at a specific can be obtained using the "pure-precession" ansatz (see also discussion in 



171]), which assumes that Si and S2 associated with the gyroscope always stay in the same 



plane as = — e^, and precess with the same angular velocity 0^6?,. More specifically, we 
can use the Euler angles in Frame I to write 

51 = -(sin 6*1 sin0ie^ — sin 6^1 cos0iej^ + cos6'ie^), (Ala) 

52 = -(sin^2 sin02e^ - sin02 cos02ej, + cos6'2e^). (Alb) 

2 ^ 

The ansatz assumes that 61 = 62 = 0, <pi = <p2 = 4>p, and 02 — 0i = tt (the last relation 
can be seen from the initial configuration at the neutrino sphere with Si = — S2 = e^/2 
corresponding to 61 = 26^ and 62 = tt — 26y). Using this ansatz along with conservation of 
5*2 and Eqs. fl20al) and fl20bl) . we obtain 

cos 9i + a cos 6^2 = (1 — a) cos 29^, (A2a) 



bp sin 01 = jj,^ sin ^1 + ^ sin (6^1 + 6^2), (A2b) 



0psin^2 = -/ivsin^2 + it sin(6'i + 62). (A2c) 

For any specific /i, the above equations can be solved to give ^i, Q21 and 0p. Then we can 
calculate the corresponding a, Q, and using Q\ and Q2- The above procedure gives a, Q, 
Qp, and 0p for the gyroscope at a specific \i that are indistinguishable from those obtained 
by the procedure discussed in Sec. IIII B[ 
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Appendix B: Initial Conditions for the Neutrino Gyroscope 



The initial conditions for the neutrino gyroscope at the neutrino sphere are 



S(0) 



I — a 



Q(0) = ii^ei + ^H.. 

Noting that Hv = — sin26've^ — cos26've^ = — e^, we obtain from the above equations 



Sz = — cos 26^ 



Q 



a 



1 + a fi^ 



cos 26, 



2 11 
1 — a f 1 + a fi 



sin 26, 



1/2 



^cos2^, 
/i 



(Bla) 
(Bib) 

(B2a) 
(B2b) 
(B2c) 



2Q V 2 

In terms of the dynamic variables 6 and 0, the initial conditions for the neutrino gyroscope 
can be chosen as 6'o = 0, 0o = 0, and 



1 fl + a ~ fx 

cos 6o = \ — - — cos 26'v 

V V 2 



50 = Atv 



1 — a 
l + a 



(B3a) 
(B3b) 



where the subscript "0" indicates the initial moment t = and Eq. flB3bl) is obtained from 
Eq. (12 lap at t = 0. For a constant /i, the precession frequency at 6 = 6p, where Ves{6) 
reaches its minimum, is given by Eq. ( l3Ti) as 



fia ± y/{fi(ry - Afifi^Q cos 6p _ J (^f^/ cos^p) - (fivQ/cr), 



2 cos 6r, 



(B4) 



where the approximate equalities apply for fi^/ fi with the upper and lower expressions 
corresponding to the plus and minus signs in front of the square root, respectively. It can 
be shown that the upper expression of (pp is unphysical as it cannot satisfy conservation of 
Sz- For the physical value of 0p ~ fi^Q/a, conservation of gives 

2(1 + a) /x^ 



6p ^ 26, 



1 + 



'1 — ay jjL 



(B5) 



The above expression of 6p is to the first order in /iv//i and 6,. To the same order, we have 



6. ^26Al 



2 /Xy 

1 + a /i 



(B6) 
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The amplitude of nutation around 6 = 6p is then 



ri = 0p-0o^ 



(l + a)(l-a)2 /i 
The number density of z/g at the neutrino sphere is 



FoT fi = 2V2GFn,^{R,), 



1.66 X 10 



IQ51 erg/s 



10 kmV MeV\ 



5.92 X 10 



-6 



fx 



5rn? 



lOMeV^ 



E 



10^^ cm ^ 



cm 



(B7) 



(B8) 



^3 X 10-3 eV^ 

Taking 0^ = 10-^ Sm"^ = 3 x IQ-^ eV^, 2^1^ /5m^ = 1/9 MeV-\ a = 2/3, and n^^R^) 
1.66 X 10^^ cm-3, we have 2.28 x 10"^. 



(B9) 



Appendix C: The Neutrino Gyroscope at the Critical Point 

The critical point at /i = separates the evolution of the neutrino gyroscope into 
two regimes: the sleeping-top regime with essentially pure precession but little nutation at 
fi ^ ficr and the other with both precession and nutation at yU < /icr- As ^ ^ 1 in the 
sleeping-top regime, we derive /icr assuming 9y = 0. Let the gyroscope start with 6* = at a 
constant fi. Its fixed parameters are 

Sz = a = (Cla) 

g = - ^, (cib) 



The motion of the gyroscope is governed by 



^gyro = + /^vQ. (Clc) 



-sin^^ = (t(1 -cos^), (C2a) 
^(^^ + 02sin^^) =/ivQ(l-cos0), (C2b) 

^ Laj 



which are obtained by rewriting Eqs. fl25cp and fl26bp . Assuming that > is allowed. 



we can find ^max, the maximum value of ^, by setting ^ = in Eq. (]C2b|) . Combining the 
resulting equation with Eq. ( ]C2al) . we obtain 



39 



It can be seen that when /icr^ > 4:fi^Q, the above equation has no solution for ^max > 0. 
Thus, the gyroscope remains in its initial vertical position {9 = 0) for fia"^ > Afi^Q. For a 
gyroscope starting ai fj, fj,^, this condition corresponds to /z > /Xcr, where 

4/iv 



(C4) 



Expanding = (1 — «) cos 26^^/2 to the leading order in 9^, we obtain the pure-precession 
solution at /i = from Eqs. (ISTl) . f l37ap . fl38ap . and ( I38bj) (setting = in the last 
equation) : 



a 



(0) 



■9, 



Qcr ~ Q 



(0) 



1 



8V2a 



4V2a 



p,cr 



'p,cr 



_V2a^~2/3 



0f 



+ 



4v^ ^'^7' 



where = (1 - a)/2, g(°) = (1 + a)/2 - fx^/^icr, and 



(C5a) 
(C5b) 
(C5c) 
(C5d) 

/icr(T^°V2- For a = 2/3 and 



9^ 



10 ^, we have 6'pcr ~ 8-24 x 10 which agrees with the numerical result very well. 



Using the above results at the critical point and Eq. fl30b|) . we obtain 

x/3 



-/icrO" 



p,cr • 



(C6) 



We note that the behavior of the precession and nutation frequencies at the critical point as 
shown in Fig. |9]is unique to the IH. In contrast, both the precession and nutation frequencies 
increase smoothly with /x for the NH. 



[1] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 070521 (2010). 
[2] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978). 

[3] S. P. Mikheyev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 

[4] H. Duan, G. M. Fuller, and Y.-Z. Qian, Annu. Rev. Nucl. Part. Sci. 60, 569 (2010), 

arXiv:1001.2799 [astro-ph]. 
[5] G. M. Fuller, R. W. Mayle, J. R. Wilson, and D. N. Schramm, Astrophys. J. 322, 795 (1987). 



40 



[6] J. Pantaleone, Phys. Lett. B 287, 128 (1992). 

[7] G. Sigl and G. G. RafFelt, Nucl. Phys. B406, 423 (1993). 

[8] A. B. Balantekin and H. Yiiksel, New J. Phys. 7, 51 (2005), astro-ph/0411159. 

[9] G. M. Fuller and Y.-Z. Qian, Phys. Rev. D 73, 023004 (2006), astro-ph/0505240. 
[10] H. Duan, G. M. Puller, and Y.-Z. Qian, Phys. Rev. D 74, 123004 (2006), astro-ph/0511275. 
[11] H. Duan, G. M. Fuller, J. Carlson, and Y.-Z. Qian, Phys. Rev. D 74, 105014 (2006), astro- 
ph/0606616. 

[12] G. G. RafFelt and A. Y. Smirnov, Phys. Rev. D 76, 081301(R) (2007), arXiv:0705.1830. 
[13] G. Fogh, E. Lisi, A. Marrone, and A. Mirizzi, JCAP 0712, 010 (2007), arXiv:0707.1998. 
[14] B. Dasgupta, A. Dighe, G. G. Raffelt, and A. Y. Smirnov, Phys. Rev. Lett. 103, 051105 

(2009), arXiv:0904.3542. 
[15] A. Mirizzi and R. Tomas, preprint (2010), arXiv: 1012. 1339 [hep-ph]. 

[16] S. Hannestad, G. G. Raffelt, G. Sigl, and Y. Y. Y. Wong, Phys. Rev. D 74, 105010 (2006), 
astro-ph/0608695. 

[17] H. Duan, G. M. Fuller, J. Carlson, and Y.-Z. Qian, Phys. Rev. D 75, 125005 (2007), astro- 
ph/0703776. 

[18] H. Duan, G. M. Fuller, and Y.-Z. Qian, J. Phys. G 36, 105003 (2009), arXiv:0808.20 [astro-ph]. 

[19] G. G. RafFelt and A. Y. Smirnov, Phys. Rev. D 76, 125008 (2007), astro-ph/0709.4641. 

[20] G. G. Raffelt, preprint (2011), arXiv:1103.2891 [hep-ph]. 

[21] S. Calais and C. Volpe, preprint (2011), arXiv: 1103.5302 [astro-ph]. 



41 



